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ABSTRACT 

Near  surface  measurements  of  turbulent  temperature  and 
velocity  fluctuations  in  the  ocean  are  obscured  by  relatively 
much  larger  ordered  fluctuations  due  to  surface  gravity  waves. 
Analysis  of  measured  velocity  signals  in  the  near  surface 
region  indicated  that  the  ratio  of  turbulent  to  wave-induced 
kinetic  energy  was  consistently  less  than  10~  .   To  a  first 
approximation,  wave^induced  fluctuations  should  not  contribute 
to  turbulent  transport.   However,  small  instrument  errors  due 
to  misalignment,  phase  response,  or  directional  response  may 
cause  significant  wave^induced  error  in  direct  flux  measure- 
ments.  The  effect  on  turbulent  flux  measurements  of  misalign- 
ment of  the  velocity  sensors  was  examined.   It  was  found  that 
misalignment  resulted  in  apparent  phase  shifts  in  the  measured 
wave-induced  signals,  causing  erroneous  contributions  to  the 
calculated  flux.   For  alignment  errors  of  3.6°,  the  error  in 
calculated  momentum  and  heat  fluxes  could  have  been  up  to 
500%  abd  110%,  respectively.   Momentum  flux  measurements  made 
in  deep  water  (with  respect  to  the  length  of  surface  waves) 
were  found  to  be  less  sensitive  to  alignment  errors.   The 
errors  in  this  case  were  functions  of  the  alignment  of  the 
flowmeter  axes  with  the  predominant  wave  direction. 
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It   INTRODUCTION 

Turbulent  processes  in  the  upper  ocean  are  not  yet  fully 
understood,   Phillips  (1966),  Kitaigorodskii  (1971)  and  others 
have  pointed  out  the  need  for  near  surface  measurements  of 
*   turbulence  and  turbulent  transfer.   Large  coordinated 
experiments  such  as  BOMEX  have  been  designed  to  measure 
turbulence  in  the  Marine  Boundary  Layer  of  the  atmosphere 
(Holland,  1972),  but  complementary  measurements  in  the  near 
surface  layer  of  the  ocean  have  been  scarce.   The  recently 
completed  GATE  should  greatly  add  to  such  data.   Several 
projects  (NORPAX,  JASIN,  GARP)  are  currently  being  planned 
to  measure  turbulence  in  what  is  generally  termed  the  mixed 
layer  of  the  ocean. 

The  dynamics  of  the  near  surface  layer  (in  the  absence 
of  tropical  storms)  are  dominated  by  surface  gravity  waves. 
The  ordered  motion  of  the  waves  is  usually  not  of  interest 
when  making  turbulence  measurements  because  they  do  not 
contribute,  at  least  to  a  first  approximation,  to  the  vertical 
transfer  of  heat  and  momentum  through  the  water  column,   As 
a  result,  near  surface  turbulence  measurements  seek  to 
resolve  what  are  essentially  second  order  quasi-random 
motions  against  a  background  of  first  order  coherent  noise. 

The  Naval  Postgraduate  School  has  conducted  a  series  of 
turbulence  measurements  using  a  fixed  array  of  sensors 
mounted  on  the  Naval  Undersea  Center  oceanographic  research 
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tower  off  San  Diego  (Hagen,  1974),   An  attempt  was  made  to 
measure  turbulent  heat  and  momentum  fluxes  from  simultaneously 
recorded  time  series  of  orthogonal  velocity  components  and 
temperature.   It  was  found  that  wave-induced  fluctuations 
may  have  caused  significant  errors  in  the  measurements.   In 
view  of  the  current  emphasis  on  near  surface  turbulence 
measurements,  an  analysis  of  wave-induced  measurement  errors 
was  considered  appropriate.   Only  through  such  analyses  can 
realistic  experiments  be  designed. 

The  analysis  of  wave-induced  errors  in  turbulence  measure- 
ment is  carried  out  first,  followed  by  a  description  of  the 
experiment  and  the  results  obtained. 

II.   WAVES,  TURBULENCE  AND  NOISE 

In  the  upper  layer  of  the  ocean,  turbulence  is  superim- 
posed on  ordered  motions  related  to  surface  gravity  waves. 
The  contribution  of  the  waves  to  a  measured  signal  may  be 
quite  large.   Data  taken  from  the  NPS  experiments  were 
examined  as  to  the  relative  magnitudes  of  wave-induced 
signals,  turbulence  and  noise.   The  data  were  presumed  repre- 
sentative of  open  ocean  as  well  as  coastal  near  surface 
environments. 

Fluctuations  of  horizontal  and  vertical  velocity  (u,  w) 
and  temperature,  0,  can  be  thought  of  as  composed  of  wave- 
induced  and  turbulent  parts: 
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u  =  U  +  u' 
w  =  W  +  wf 

e  =  e  +  e*  (i) 

where  upper  case  letters  stand  for  wave-induced  fluctuations 
and  primed  letters  for  turbulent  fluctuations.   With  certain 
assumptions  it  is  possible  to  separate  the  wave-induced  and 
turbulent  parts  of  equation  C.1)  and  estimate  their  relative 
magnitudes.   For  moderate  wave  conditions,  it  is  reasonable 
to  assume  that  the  wave-induced  spectral  velocity  components 
can  be  calculated  from  the  wave  spectrum  using  linear  wave 
theory.   The  success  with  which  wave  velocities  are 
described  by  linear  wave  theory  can  be  gauged  by  applying 
the  appropriate  linear  transfer  function  to  the  surface 
displacement  spectrum  and  studying  how  well  the  calculated 
and  measured  velocity  spectra  compare.   Figure  (1)  shows  the 
measured  vertical  velocity  spectrum  at  6.6  meters  depth  and 
that  calculated  from  the  surface  displacement  spectrum  using 
linear  wave  theory.   Figure  2  shows  the  coherence  and  phase 
functions  obtained  from  the  cross  spectrum  between  measured 
wave  and  velocity  records.   The  maximum  coherency  squared 
(hereafter  referred  to  simply  as  the  coherence)  occurs  at  a 
period  of  5.9  seconds  and  is  equal  to  0.9948.   The  phase 
function  shows  the  waves  leading  the  vertical  velocities  by 
approximately  90°  throughout  the  surface  wave  band.   Figures 
1  and  2  clearly  demonstrate  the  applicability  of  linear  wave 
theory  to  wave-induced  motions. 
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Figure  1,   Comparison  between  measured  vertical  velocity 
spectrum  and  the  vertical  velocity  spectrum 
calculated  from  the  waves  using  linear  theory 
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Figure  2,   Coherence  and  phase  functions  calculated  from  the 
cross  spectrum  between  surface  displacement  and 
vertical  velocity. 
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It  is  assumed  that  wave-induced  and  turbulent  motions 
are  independent  (Phillips,  1966)  (Seitz,  1971),   The  inter- 
actions between  wave  motions  and  turbulence  contribute  to 
the  motion  only  at  the  third  order  (Kinsman,  1965),   There- 
fore, the  independence  assumption  is  justified  except  for 
very  steep  waves. 

Using  the  two  preceding  assumptions,  the  measured  velocity 
spectrum  can  be  separated  into  wave-induced  and  turbulent 
parts  using  the  experimentally  determined  coherence  between 
the  wave  displacement  and  velocity  signals  (Benilov  and 
Filyushkin,  1970,  see  Appendix  A),   The  turbulence  spectrum 
is  then  given  by 

S  ,(f)  =  S  (f)Cl  -  Y2   )  (2) 

w*   '    w        '  nw  v  ' 

and  the  wave-induced  spectrum  by 

Sw(f)  =  Sw(f)  -  Sw,(f),  (3) 

The  variance  of  velocity  due  to  wave-induced  or  turbulent 
motions  is  obtained  by  summing  their  respective  spectral 
components.   Table  1  lists  experimentally  determined  velocity 
variances  from  the  NPS  data,   The  horizontal  wave  velocities 
are  larger  than  the  vertical  velocities  due  to  the  shallow 
water  environment.   The  magnitude  of  horizontal  "turbulent" 
variance  may  not  all  be  attributable  to  turbulence  due  to 
the  effect  of  directionality  of  horizontal  wave  motions  on 
the  coherence  (Yefimov  and  Khristoforov,  1971).   Table  1 
indicates  that  the  intensity  of  turbulent  motions  decreases 
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RUN 


DEPTH 


VELOCITY  VARIANCE  (cm2sec"2) 


8,04 

-2 

u 

-2 

w 

17 

176.6 

28,5 

18 

12.54 

228,7 

14,0 

19 

9.44 

229.6 

35.6 

22 

6.62 

310,9 

157.7 

26 

5,84 

577.9 

306 , 2 

28 

10.95 

161,1 

11,8 

28a 

10.95 

133.4 

10.8 

U2 

W2 

u'2 

w'2 

168,4 

25.5 

7.7 

2.9 

208.2 

12.4 

20,2 

1.6 

210.5 

33.1 

18.6 

2.4 

296,0 

150,8 

14.2 

6,5 

555,8 

294.1 

20.7 

11,3 

152,2 

11.1 

8.5 

0.7 

126.0 

10.1 

7.1 

0.7 

TABLE  I ,   Velocity  variance  calculated  from  wave-induced  and 
turbulent  velocity  spectra. 


RUN 

DEPTH 

TEMPERATURE  VARIANCE 

(°c2) 

7.98 

e2 

02 

e2 

17 

0,0787 

0.0115 

0.0672 

18 

12,48 

0.0552 

0,0026 

0,0526 

19 

9,38 

0,0089 

0,0016 

0.0073 

22 

6,56 

0,0192 

0,0067 

0,0125 

26 

5.78 

0.0178 

0,0051 

0.0127 

28 

10.89 

0.0126 

0.0019 

0.0107 

28a 

10.89 

0.0485 

0.0095 

0.0390 

TABLE  II.   Temperature  variance  calculated  from  wave-induced 
and  turbulent  temperature  spectra. 
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with  depth  roughly  as  the  intensity  of  the  wave  motions.   The 
ratio  of  turbulent  to  wave^induced  energy  is  seen  to  be 
consistently  less  than  10^  . 

An  inseparable  part  of  the  turbulent  variance  is  that  due 
to  instrument  noise.   The  advertised  broadband  noise  level 

for  the  flowmeters  used  CMarsh-McBirney  Model  711)  was  1.0 

2   -2 

cm  sec   ,   The  noise  level  was  checked  by  recording  the 

instrument  output  while  the  flowmeters  were  at  rest  in  still 
water  and  was  visually  found  to  agree  with  that  advertised. 
Additional  noise  can  be  introduced  in  the  form  of  turbulence 
created  by  the  obstruction  of  the  flow  field  by  the  sensor. 
At  the  small  velocities  measured  this  effect  can  probably  be 
ignored  (NOIC,  1974). 

Making  the  additional  assumption  that  temperature  can  be 
treated  as  a  passive  scalar,  the  assumptions  made  with  regard 
to  velocity  fluctuations  should  apply  equally  well  to  the 
wave-induced  and  turbulent  temperature  fluctuations  of 
equation  (1).   The  wave-induced  temperature  fluctuations  can 
then  be  written  in  terms  of  the  vertical  displacement  of 
some  local  temperature  gradient  relative  to  the  point  of 
measurement  (Thornton  and  Boston,  1974): 


e  =  zi|S=  *  »f°  ,  (4) 


where  Z  is  the  wave-induced  vertical  displacement  and 
39b/ 8z  is  the  temperature  gradient.   Equation  (4)  implies 
that  the  coherence  between  0  and  the  surface  displacement 
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is  1,0  and  that  the  method  of  Appendix  A  can  be  used  to 
separate  wave-induced  and  turbulent  temperature  components. 
Table  2  shows  the  wave-induced  and  turbulent  temperature 
variances  calculated  from  the  NPS  data, 

III,   APPLICABILITY  OF  TAYLORS  HYPOTHESIS 

An  associated  problem  of  turbulence  measurements   in 
the  near  surface  layer,  although  not  directly  tied  to  flux 
measurements,  deals  with  the  applicability  of  Taylor's 
hypothesis.   Experimental  turbulence  data  is  usually  in  the 
form  of  a  time  history  of  turbulent  fluctuations  measured 
at  a  point  or  with  a  towed  sensor.   The  spectra  calculated 
from  such  measurements  are  then  defined  in  frequency  space. 
Unfortunately  the  theoretical  treatment  of  turbulence  is 
usually  confined  to  wave  number   space,  so  that  a  trans- 
formation is  employed.   Turbulent  frequencies  are  not  uniquely 
related  to  their  wavenumber,  as  is  the  case  for  wave  motions 
for  which  a  dispersion  relationship  exists.   Instead,  use 
is  made  of  Taylor's  hypothesis,  or  the  frozen  turbulence 
approximation.   Taylor's  hypothesis  states  that  if  a  sensor 
is  moved  rapidly  enough  through  a  turbulent  field,  the  time 
dependent  changes  in  the  field  are  negligible,  with  the 
result  that  it  is  the  spatial  variations  that  are  being 
measured,   A  transformation  from  time  to  space  can  therefore 
be  made  by  substituting  x/U  for  t,  where  U  is  the  velocity  of 
the  sensor  (Tennekes  and  Lumley,  1972),   For  fixed  point 
measurements  the  transformation  can  still  be  made,  providing 
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that  the  turbulent  field  is  being  swept  by  the  sensor  fast 
enough,  as  by  a  mean  wind  or  current.   The  criterion  for  the 
applicability  of  Taylor's  hypothesis  is  that  the  relative 
velocity  of  the  sensor  with  respect  to  the  fluid  be  much 
larger  than  the  magnitude  of  turbulent  velocity  fluctuations: 

u7  U  <<  1, 

Meteorologist  never  seem  to  have  trouble  making  the 
transformation  because  mean  wind  speeds  are  on  the  order  of 
meters  per  second  while  small  scale  turbulent  fluctuations 
are  on  the  order  of  centimeters  per  second.   Mean  flow 
velocities  in  the  ocean  are  usually  measured  in  centimeters 
per  second,  as  are  turbulent  fluctuations.   In  the  ocean  it 
is  necessary  to  verify  whether  or  not  Taylor's  hypothesis 
is  valid,  or  to  tow  the  sensor  through  the  water  at  a  speed 
considered  sufficient  to  satisfy  the  approximation.   The 
latter  method  has  been  used  most  often  to  date. 

Seitz  (1971)  and  Bowden  (1963)  concluded  Taylor's 
hypothesis  was  valid  in  a  tidal  estuary  for  frequencies  down 
to  0,25  Hz  when  the  mean  convective  velocity  was  on  the  order 
of  25  cm/sec'  .   The  NPS  experiment  did  not  include  a  good 
measurement  of  mean  velocity,  but  turbulent  fluctuations  were 
on  the  order  of  0.3  cm  sec'"  per  spectral  estimate  (Af  = 
0.0067  Hz).   A  mean  velocity  of  only  3.0  cm  sec"   could 
conceivably  be  sufficient  to  satisfy  Taylor's  hypothesis  in 
this  case. 
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High  mean  velocities  are  seldom  encountered  in  the  open 
ocean.   Even  where  surface  currents  may  be  large  enough,  the 
measuring  platform  usually  moves  with  the  current  so  that 
relative  velocities  might  be  negligible.   In  such  cases,  it 
is  doubtful  whether  Taylor *s  hypothesis  can  be  applied. 

IV,   THEORY  OF  TURBULENT  FLUX  MEASUREMENT 

If  the  Navier-Stokes  equation  is  rewritten  using  mean 
and  fluctuating  velocities  and  then  time  averaged,  one  of 
the  resulting  terms  is  the  Reynold's  stress  tensor  (Phillips, 
1966): 


-  p  UjU     i,  j  =  1,  2,  3  (5) 

■J 

were  u.  are  velocity  fluctuations  and  the  overbar  denotes  a 

time  average.   The  off  diagonal  terms  {if j ) ,  represent  a 

turbulent  transport  of  momentum  and  are  covariances  of 

orthogonal  velocity  fluctuation  components.   If  the  motion 

can  be  considered  horizontally  homogeneous,  the  downstream 

component  of  vertical  momentum  flux  is  given,  in  conventional 

notation,  by 

T 

u(t)w(t)  dt,  (6) 

An  analogous  expression  for  vertical  heat  flux  can  be 

derived  from  thermodynamic  relations  and  is 

T 

9(t)w(t)  dt,  (7) 
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where  G  represents  a  fluctuating  temperature,   The  covariance 
between  two  fluctuating  quantities  can  also  be  expressed 
(Bendat  and  Piersol,  1971) 

uw  =  R   CO)  =  2     Co   (f)  df  ,  (8) 

uwv  '      /     uw  v  ' 

J  o 

where  Co   (f)  is  the  real  (_in   phase)  part  of  the  u-w  cross- 
spectrum,   The  spectral  representation  is  the  most  popular 
because  turbulence  theory  has  envolved  in  the  spectral  domain 
Experimental  determination  of  heat  and  momentum  flux  can  be 
made  by  recording  fluctuating  velocity  components  u  and  w, 
and  temperature  fluctuations,  9,  at  a  point  and  calculating 
the  covariance  or  cospectrum  as  above. 

Alternate  methods  for  estimating  turbulent  heat  flux 
are  the  dissipation  method  used  by  Williams  and  Gibson 
(1974)   and  the  profile  method.   The  latter  is  essentially 
an  integration  over  time  of  an  assumed  equation  of  conser- 
vation for  heat  (Voskanyan,  et ,  al , ,  1969)(Kaiser  and 
Williams,  1974).   Both  these  methods  have  the  advantage  that 
only  temperature  measurements  are  necessary,  while  the  eddy 
correlation  method  discussed  above  requires  simultaneous 
temperature  and  velocity  measurements.   An  advantage  of  the 
eddy  correlation  method  is  that  the  kinematic  regime 
associated  with  turbulent  heat  flux  can  be  analyzed. 


Measures  only  the  eddy  diffusion  coefficient  for  heat 
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V,   WAVE  EFFECTS  ON  FLUX  MEASUREMENT 

The  wave-induced  contribution  to  velocity  signals  in  the 
near  surface  region  has  been  shown  to  be  much  larger  than 
the  turbulent  contribution.   The  effect  of  wave  motions  on 
the  measurement  of  turbulent  flux  is  analyzed  below. 

The  covariance  between  u  and  w  in  a  region  of  wave 
influence,  assuming  turbulent  and  wave^induced  motions  are 
statistically  independent,  can  be  written: 


uw  =  UW  +  u'w' .  (9) 

For  simplicity,  u  refers  to  only  one  horizontal  velocity 
component.   The  first  term  is  the  contribution  of  the  wave- 
induced  motions  and  the  second  the  contribution  due  to 
correlated  turbulent  velocity  fluctuations.   A  similar 
expansion  into  wave-induced  and  turbulent  parts  can  be  made 
for  temperature  and  vertical  velocity  covariance: 


6w  =  0W  +  9*w* ,  (10) 

Since  it  has  been  shown  that  linear  wave  theory  closely 
approximates  observed  data,  U  and  W,  as  well  as  0  and  W 
(equation  4)  are  in  quadrature  and  their  contributions  to 
the  respective  covariances  should  be  zero.   In  terms  of  the 
cross  spectrum: 


S   (f)  =  Co   (f)  -  iQ   (f),  (11) 

uwv  '  uwv  '     niwv  '  ' 


wave  motions  should  contribute  only  to  the  imaginary  part, 
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Recalling  the  definition  C8),  equation  (9)  can  be  written 

uw  =  Re  [2/   Suw(f)  df]=  Re  [uhv'  ~  i(U  W  +  6)]  ,    (12) 
-'  o 

The  contribution  of  the  wave  motions  appears  only  in  the 
(imaginary)  quadrature  spectrum  because  U  and  W  are  90°  out 
of  phase.   The  quadrature  contribution  of  the  turbulence  (6) 
will  henceforth  be  neglected  compared  to  UW,   The  simple 
overbar  denotes  RMS  units,   If,  for  some  reason,  the  measured 
phase  between  U  and  W  should  differ  from  90°  a  component  of 
one  signal  will  appear  in  phase  with  the  other  and  contribute 
to  the  cospectrum,   Such  a  contribution  results  in  an  error 
in  the  calculated  flux.   This  error  might  be  ignored  for 
small  phase  differences  were  it  not  for  the  fact  that  the 
wave-induced  signal  is  so  much  larger  than  the  turbulent 
signal. 

Real  phase  shifts  between  U  and  W  due  to  non-linear 
effects  have  been  postulated.   For  example  Yefimov  and 
Khristoforov  (1972)  included  viscosity  effects  in  the  equations 
of  motion  to  show  that  a  phase  shift  occurs  at  second  order 
for  steep  waves.   The  phase  shift  attenuates  rapidly  with 
depth  and  increases  with  increasing  frequency.   Consistent 
frequency  dependent  phase  shifts  were  conspicuously  absent 
from  even  our  shallowest  velocity  records  in  the  entire  wind 
wave  range  of  frequencies  (see  for  example  figure  1). 

Other  than  unknown  higher  order  or  weak  nonlinear  effects, 
the  only  cause  for  the  phase  between  horizontal  and  vertical 
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wave-induced  velocities  to  differ  from  90°  is  measurement 
error  in  directional  response,  frequency  response,  or  improper 
instrument  orientation,   Non  ideal  directional  response 
will  be  examined  first, 

A.   FLOWMETER  DIRECTIONAL  RESPONSE 

Expressing  the  true  harmonic  velocity  component  as 

V  =  | V |   cos  o  t  (12) 

the  ideal  response  of  a  vertically  oriented  flowmeter  is 
pictured  in  polar  and  cartesian  form  in  figure  3a.  If  the 
flowmeter  response  deviates  from  the  cosine  law,  the  measured 
signal  will  be  in  error  by  some  amount  e(at)  which  depends 
on  the  magnitude  and  phase  (i,e,  direction)  of  V.   The 
measured  signal  becomes 

V*  =  | V |   cos  a t  +  | V |  e(at).  (13) 

z^ot)   must  be  periodic  in  at  and  can  thus  be  decomposed  into 

Fourier  components  which  are  all  harmonics  of  a  but  are  not 

necessarily  in  phase  with  V, 

m 
e(at)  =  Y.     z      cos(nat  +  <(>).  (14) 

n=l   n  n 

There  are  three  forms  in  which  the  error  can  appear. 

First,  <b      can  be    zero  by  virtue  of  the  fact  that  the  flow- 
rn 

meter  will  not  respond  to  motions  perpendicular  to  its  axis. 
The  poor  directional  response  of  the  flowmeter  in  this  case 
has  the  effect  of  broadening  the  bandwidth  of  the  measured 
velocity  signal  by  inclusion  of  the  harmonics  of  e(at). 
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For  flowmeters  with  noticeably  poor  directional  response  (as 

in  figure  3b),  this  may  create  difficulties  in  analyzing 

velocity  power  spectra  for  slopes  indicative  of  an  inertial 

subrange  above  the  wave  frequency  band.   The  electromagnetic 

flowmeters  used  here,  and  others  tested  by  NOIC  (1974),  showed 

a  typical  first  harmonic  error  response  on  the  order  of  2  to 

4  percent  of  the  ideal  response  amplitude.   Table  I  indicates 

that  the  magnitude  of  turbulent  fluctuations  was  on  the  order 

of  10  percent  of  the  wave  induced  amplitudes.   Contamination 

of  the  high  frequency  turbulent  signal  by  harmonics  of 

e(at)  might  therefore  be  possible. 

As  long  as  the  flowmeter  does  not  respond  to  perpendicular 

flow,  however,  the  measured  signal  must  be  in  phase  with  the 

actual  motion.   If,  for  example,  the  covariance  is  calculated 

between  W*  and  another  signal  with  which  W  should  be  in 

quadrature  (say  U),  the  two  measured  signals,  including  all 

harmonics  of  e(at),  are  still  in  quadrature.   By  virtue  of 

their  orthogonality  there  is  no  erroneous  flux  contribution. 

Surface  waves  usually  occupy  a  frequency  band  of  at  least 

two  octaves  and  there  is  a  possibility  that  the  first  harmonic 

of  the  error  signal  in  W*  is  not  in  quadrature  with  the  U 

signal  at  that  frequency.   Since  the  phase  between  spectral 

wave  components  is  random,  however,  the  average  contribution 

can  be  expected  to  be  zero. 

The  second  case  of  (14)  is  when  the  <J>   are  zero  but 

n 

e(+  ^)  is  not  because  of  a  D.C.  bias  in  the  system.   Errors 
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of  this  sort  are  inconsequential  because  DtC,  values  are 
normally  subtracted  out  of  the  fluctuating  signal  before 
calculations  are  made. 

The  last  schematic  in  figure  3  is  for  the  case  of  (J)  ^  0. 
This  kind  of  response  occurs  when  the  measured  signal  at 

at  =  +  ■«  differs  from  that  at  at  =  r  ^,   In  other  words,  the 

meter  response  to  a  velocity  component  perpendicular  to  its 
axis  depends  on  the  direction  of  that  velocity  component. 
An  error  of  this  type  does  result  in  a  phase  shift  of  the 
measured  signal  with  a  consequent  wave-r-induced  contribution 
to  the  calculated  covariance.   Ducted  type  flowmeters  in 
particular  can  respond  erratically  at  very  large  angles  of 
attack.   Smith  (1972)  circumvented  this  difficulty  by 
tilting  his  ducted  impeller  current  meter  array  from  the 
principal  gravitational  axis  while  measuring  turbulent  flux 
under  Arctic  ice.   Unfortunately  in  the  near  surface  region 
of  the  open  ocean,  the  three  dimensional  wave  motions  have 
no  preferred  direction  and  this  method  does  not  work. 

B.   FLOWMETER  FREQUENCY  RESPONSE 

There  is  an  inherent  phase  lag  in  the  response  of  any 
real  sensor.   Instruments  for  measuring  turbulence  are  usually 
chosen  to  exhibit  a  flat  amplitude  response  over  the 
frequency  band  being  measured.   The  phase  response  is  seldom 
reported.   Measurements  using  a  variety  of  instruments, 
whose  outputs  are  subsequently  compared,  are  particularly 
susceptible  to  phase  errors  due  to  synchronization  problems 
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and  non-uniform  phase  response  among  the  instruments.   Later 
results  will  point  out  how-  extremely  sensitive  correlation 
measurements  in  the  layer  of  wave  influence  are  to  miniscule 
errors  in  the  measured  phase.   Care  must  be  taken  to  ensure 
that  all  instruments  have  extremely  short  time  constants  or 
at  least  that  the  phase  response  of  all  instruments  is  the 
same. 

C.   FLOWMETER  MISALIGNMENT 

The  other  source  of  phase  error  in  the  measured  velocity 
signals  to  be  examined  is  that  due  to  improper  orientation 
of  the  velocity  sensors  with  respect  to  the  natural  axes  of 
wave  motion.   Since  horizontal  wave  motions  are  directional, 
the  orientation  of  the  natural  horizontal  axes  cannot  be 
specified  except  that  they  are  orthogonal  to  the  local 
gravity  vector.   For  each  small  band  of  the  wave  spectrum, 
however,  it  can  be  assumed  that  the  direction  is  constant 
for  a  reasonable  length  of  time.   With  this  in  mind, the  three 
dimensional  field  of  wave  motions  at  a  point  can  be  represented 
by  a  Fourier  sum  of  two  dimensional  motions,  namely: 

m 
W  -  I  lwnl  cos(ant  +  <f>n) 

m 
U  =  £  |U  |sin(a  t  +  «j>  )  cos6n  (15) 

n=l 

m 

V  =     7    |U    |sin(a   t   +   $    )    sin3 
*--  '    n1         v    n  Yn  n 

n=l 
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where  3  are  the  angles  between  the  horizontal  wave  component 

directions  and  some  arbitrary  reference  direction. 

Choosing  one  Fourier  component  and  setting  3  equal  to 

zero  allows  the  measurement  of  the  motion  to  be  simply 

analyzed  by  assuming  that  the  flowmeter  axes  are  in  the  plane 

of  motion  defined  by  U   and  W  .   Figure  4  shows  the  axes  of 

n      n 

the  wave  motion  and  two  flowmeter  component  axes  which  are 
tilted  by  an  arbitrary  angle.  From  figure  4,  the  measured 
signals  U*  and  W*  are  given  by 


U*  =  |U   sinat  cosa^  +  |W   cosat  cosa01 
n    '  n  I  11    '  n '  31 


W*  =  |U   sinat  cosaiC  +  W   cosat  cosa00 
n    l  n  i  13    '  n  '  33 


(16) 


The  effect  of  misalignment  of  the  velocity  sensors  on 

turbulent  flux  measurements  can  now  be  calculated  by  taking 

the  covariance  between  U*  and  W*  and  0  .   Recognizing  that 

n      n      n       &     & 

Fourier  wave  components  are  independent,  the  covariances 


U*W*  and  0  W*  can  be  summed  over  all  wave  components  to  give 
n  n      n  n  *  & 

the  total  covariance  contribution  for  the  entire  wave 

spectrum  (still  assuming  all  3   =0).   In  terms  of  RMS 
amplitude  this  results  in 

tt2  .        =2 


U*w*  =  U   cosan  sina33  -  W   cosa33  sina1:,         (17) 


0W*  =  0U  sina33  (18) 

Equations  (17)  and  (18)  are  the  errors  due  to  wave  motion  in 
measured  turbulent  flux  if  the  flowmeter  axes  are  misaligned 
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Figure  4.   Measurement  of  two  dimensional  wave  motion  by 

flowmeter  in  the  plane  of  motion.   If  the  flowmeter 


axes  are  orthogonal : 
a13  =  90  ~  a33' 


a 


11 


=  a 
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a31  =  9°  +  a33' 
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From  equations  (9)  and  (.10),  the  percent  error  introduced 
into  the  flux  measurement  is 


tt2 


u*w*   -  u'w*        U*W*        U   cos   1:Lsin33-W  cos   33sin   11 


6w*    -   9'w'        0U      since-., 

— =  dJ  (20) 

e»w»  e'w' 

The  true  value  of  the  turbulent  covariance  is  not  known,  but 
the  turbulence  correlation  coefficient  (rf)  can  be  defined: 


r*   =  u'w7(uf  wv  )  <  1  (21) 

uw        '  v   y 


r9w  =  e*wVC9*  wl  )  <  1  (22) 

Dividing  numerator  and  denominator  of  (19)  and  (20)  by  u'  w' 
or  9 '  w'   as  appropriate  gives  the  percent  error  in  terms  of 
the  tilt  angles  and  the  turbulence  correlation  coefficient. 
Equation  (20),  for  example  becomes 

9  U     sina33 

Percent  error  =  _  _ —  .  (23) 

6*w»     r' 
0w 

Equation  (23)  states  that  the  percent  error  is  inversely 
proportional  to  the  "turbulent  signal*'  to  "wave-induced  noise" 
ratio  and  the  correlation  between  turbulent  components.   It 
is  directly  proportional  to  the  tilt  angle.   Figure  5  shows 
the  percent  error  (multiplied  by  the  signal  to  noise  ratio) 
for  various  values  of  the  tilt  angle  and  turbulence 
correlation  coefficient. 
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Equation  (.19)  can  be  simplified  for  illustrative  purposes 
by  assuming  that  the  flowmeter  axes  are  orthogonal  so  that 
ot11  =  a33«   Further  assuming  that  the  tilt  angles  are  small, 
such  that  cosa11  22  l,  equation  (19)  becomes 


f2   zz2 


U*W*    CU   -  W  )  sina 
=  —  (24) 


and  the  expression  for  percent  error  can  be  written 

(U2  r.  f2)    sina33 

Percent  error  =  — — (25) 

u'w*       r' 
uw 

Figure  5  is  again  applicable  if  the  turbulent  signal  to 

wave-induced  noise  ratio  is  interpreted  as 


U2  -  W2 


Equation  (25)  shows  that,  for  the  restricted  case  described 
above,  the  error  in  deep  water  (where  W  =  U)  should  be  zero. 
An  example  of  the  momentum  flux  error  can  be  calculated 
using  the  measured  intensities  of  turbulent  and  wave-induced 
motions  from  Table  1.   Assume  a  tilt  angle  of  5°  and  a 

turbulence  correlation  coefficient,  r'  ,  of  0.2,  which  has 

'   uw'        ' 

been  borrowed  from  Bowden  and  Howe's  (1963)  results  in  an 
estuary.   The  turbulent  signal  to  wave-induced  noise  ratio 
for  run  22  is  obtained  from  Table  1  and  is  equal  to  0.064. 
Entering  figure  5,  the  product  of  the  percent  error  and 
signal  to  noise  ratio  is  44.   Solving  for  the  percent  error 
gives  a  value  of  688%, 
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Since  there  are  no  representative  values  available  for 
the  correlation  coefficient  between  turbulent  temperature  and 
vertical  velocity  fluctuations,  estimates  of  the  heat  flux 
error  rest  on  the  assumption  that  r*   is  approximately  equal 
to  r*  .   Using  the  previous  values  of  0.2  for  r*  ,  5°  for  the 

UW  or-  Qw> 

tilt  angle,  and  the  signal  to  noise  ratio  from  Tables  1  and 
2  CO, 29)  results  in  a  percent  error  for  the  heat  flux  measure- 
ment of  155%,   The  heat  flux  measurement  is  less  sensitive 
to  tilt  errors  because  of  the  larger  signal  to  noise  ratio. 

The  accuracy  with  which  flowmeters  must  be  aligned  to 
achieve  an  acceptable  percent  error  can  also  be  calculated. 
For  a  100%  error,  a  signal  to  noise  ratio  of  0,064  and  an 
assumed  turbulence  correlation  coefficient  of  0.2  the 
required  alignment  accuracy  for  momentum  flux  measurements 
is  +  0,73°,   For  the  heat  flux  measurement,  using  the  same 
values  as  the  previous  example,  the  required  alignment 
accuracy  is  3.3°, 

For  the  general  case  where  the  horizontal  axis  of  the 
flowmeter  does  not  lie  in  the  plane  of  all  orbital  motions 
(due  to  the  directionality  of  the  waves),  the  above  develop- 
ment becomes  more  complicated.   In  this  case,  let  the 
reference  direction  be  the  azimuth  of  the  horizontal  flowmeter 

axis  so  that  8   are  the  angles  the  horizontal  wave  motions 
n 

make  with  the  flowmeter  axis,  as  in  figure  6.   With  the 
coordinates  defined  as  above  all  wave  motions  are  in  the 
1-3  plane  and  V  is  zero  for  each  individual  spectral  component 
The  response  of  the  flowmeter  to  a  typical  wave  field  of 
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Figure  6.   Measurement  of  Directional  Two  Dimensional 
Wave  Motion  by  a  Flowmeter  not  in  the  Plane 
of  motion. 


35 


random  direction  is 


cosCa11)n  cos(a12)n  cos(a13), 
V*  W*)  =  CUn  V   |  (28) 

cos(a31)n  cos(a32)n  cosCc^)^ 

where  the  last  term  on  the  right  consists  of  the  directional 
cosines  between  the  unprimed  and  primed  axes.   The  angles 

(a. .)   now  include  the  effect  of  wave  direction  and  may  be 

ij  n  J 

different  for  different  Fourier  components.   Considering  a 

single  wave  component  making  an  angle  3  with  the  azimuth  of 

the  horizontal  flowmeter  axis,  the  covariances  between  U* ,  V* 

'  n   n 

and  W*  and  between  0   and  W*  are  easily  calculated 
n  n      n 

,2  T1T2 


U*w*  =  U   cosa--cosa-0  *•  W  cosa01cosa00  (29) 

n  n    n     11    13    n     31    33 


V*w*  =  -U   cosa10cosa10  +  W„  cosaQOcosaOQ  (30) 

n  n     n     12    13    n     32    33 


0  W*  =  0  U   cosaOQ 
n  n    n  n     33 


The  summation  of  wave  components  can  not  be  accomplished, 
however,  without  knowledge  of  the  direction  spectrum  of  the 
waves.   Even  then  it  is  difficult. 

Again  assuming  that  the  flowmeter  axes  are  orthogonal) 
some  statements  can  be  made  about  the  wave-induced  error. 
It  is  clear  that  the  error  contribution  to  the  U*  signal  is 

Wn  sin(a33)n 

regardless  of  the  direction  of  the  waves.   On  the  other  hand, 

the  error  contribution  to  W*  is  a  maximum  at  $^  =  0  (i.e. 

n  n 
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when  UQ  is  largest)  and  decreases  as  the  quantity  cos3  , 
The  previous  result  for  deep  water  thus  does  not  hold  and 
some  error  is  to  be  expected  in  deep  water  flux  measurements. 
The  error  for  $   =  0  in  deep  water  is  zero,  as  in  equation 
(24),   But  the  error  for  g   =  90°  is  a  maximum  given  by 

"  Wn  cosCa33)nSin(a33)n»  (31) 

since  the  U*  signal  consists  only  of  the  W  projection  onto 
the  1*  axis. 


VI,   DESCRIPTION  OF  THE  EXPERIMENT 

The  Naval  Postgraduate  School  CNPS)  has  conducted  a 
series  of  turbulence  measurements  using  a  fixed  array  of 
sensors  mounted  on  the  NUC  oceanographic  research  tower  off 
San  Diego.   The  array  was  designed  to  simultaneously  measure 
surface  elevation   temperature,  water  particle  velocity, 
conductivity  and  sound  transmission  parameters  in  a  two  cubic 
meter  volume  of  ocean.   A  detailed  description  of  the  sensors 
and  ancillary  electronics  was  given  by  Hagen  (1974),   The 
data  reported  here  consists  of  the  output  of  a  resistance 
wire  wave  gauge,  a  two-component  electromagnetic  flowmeter, 
and  a  thermistor  adjacent  to  the  flowmeter,   The  instruments 
were  all  mounted  in  a  vertical  line.   The  vertical  separation 
between  thermistor  and  flowmeter  was  six  centimeters.   Analog 
sensor  outputs  were  recorded  on  FM  multichannel  tape 
recorders  and  subsequently  passed  through  matched  analog 
lowpass  filters  before  being  digitized.   Prior  to  further 
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analysis,  all  records  were  highpass  filtered  using  a  running 
average  filter  with  a  0.00675  Hz  cutoff  frequency  to  subtract 
out  the  low  frequency  fluctuations  due  to  internal  waves. 

Data  were  obtained  in  water  approximately  twenty  meters 
deep.   A  persistent  thermocline  was  observed  at  mid-depth, 
the  mean  temperature  profile  being  characterized  by  two 
nearly  isothermal  layers  separated  by  a  relatively  sharp 
gradient.   Figure  7  shows  the  position  of  the  sensors  in 
the  water  column  and  relative  to  the  thermal  structure  as 
sampled  by  bathythermograph  (Hagen,  1974).   The  thermocline 
was  found  to  correspond  closely  with  the  12°  isotherm  in 
most  cases.   Internal  wave  activity  was  observed  during  each 
run. 

The  thermocline,  when  internal  waves  are  present,  is  a 
region  of  relatively  high  velocity  shear  combined  with  the 
large  temperature  gradient.   Velocity  shear  due  to  internal 
waves  could  not  be  determined  with  any  certainty  because  the 
internal  wave  velocities  were  only  slightly  higher  than  the 
noise  level  of  the  flowmeters.   The  maximum  apparent  shear 
was  on  the  order  of  1  cm  sec"  m~  ,   The  uncertainties  in 
low  velocity  water  motion  measurements  precluded  estimation 
of  the  Richardson  number. 

Turbulent  heat  fluxes  were  calculated  for  six  runs  made 
at  various  depths.   The  array  was  positioned  in  each  case 
such  that  the  sensors  were  near  and,  in  some  cases,  in  the 
thermocline.   As  internal  waves  passed  the  array,  the 
Lagrangian  flow  field  was  sampled  to  a  vertical  extent  of 
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about  2rv3  meters,   Thus,  the  averages  calculated  from  the 
experimental  records  were  spatial  as  well  as  time  averages. 

The  device  of  filtering  out  low  frequency  temperature 
and  velocity  data  merely  prevented  internal  wave  components 
from  biasing  the  measurements.   Since  available  averaging 
times  were  insufficient  to  provide  proper  estimates  at 
internal  wave  frequencies,  it  was  decided  to  ignore  them. 
Internal  wave  effects,  however,  were  still  observable  in  the 
filtered  records  because  of  the  above  mentioned  sampling 
of  the  Lagrangian  temperature  and  velocity  fields. 

The  calculated  turbulent  heat  fluxes  were  somewhat 
higher  than  expected.   Advection  of  heat  away  from  the 
experiment  site  was  discounted  as  a  reason  for  this  result 
due  to  the  low  current  velocities  in  the  area.   A  counter- 
balancing flux  contribution  at  the  internal  wave  frequencies 
which  had  been  filtered  out  was  considered.   The  flux  was 
calculated  from  unfiltered  records  but  nothing  conclusive 
could  be  observed  due  to  lack  of  sufficient  record  length. 
The  imposed  vertical  scale  of  twenty  meters  and  the  stable 
density  gradient  were  presumed  to  limit  turbulence  to 
frequencies  higher  than  that  of  the  internal  waves 
(Kitaigorodskii,  1973).   Attention  was  subsequently  focused 
on  the  extremely  large  effect  of  the  surface  waves  on  the 
temperature  and  velocity  records. 
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VII,   MEASURED  HEAT  FLUX 

The  heat  flux  measurements  made  from  the  NPS  array  were 
suspected  of  containing  significant  errors  due  to  flowmeter 
misalignment.   In  light  of  the  preceding  error  analysis  it 
was  considered  possible  to  make  certain  corrections  to  the 
flux  measurements  if  the  tilt  angle  of  the  vertical  velocity 
sensor  were  known.   The  cross  spectrum  between  surface  dis- 
placement and  vertical  velocity  was  available,  and  with  the 
appropriate  approximations,  a  relationship  between  the 
measured  phase  angle  and  the  actual  tilt  angle  was  derived. 

Several  fortunate  circumstances  combined  to  allow  the 
tilt  angle  to  be  calculated.   The  physical  mounting  of  the 
flowmeters  was  such  that  rotation  was  possible  only  in  the 
l'-3'  plane  Oeferring  to  figure  6).   In  addition,  coastal 
environments  are  usually  characterized  by  relatively  narrow 
banded  wave  direction  spectra  so  the  assumption  could  be  made 
that  the  wave  motions  were  essentially  in  one  direction. 
Since  the  predominant  wave  direction  was  consistently  within 
30°  of  the  azimuth  of  the  plane  of  the  l*-3'  flowmeter 
components,  small  angle  approximations  could  be  used.   Refer- 
ring again  to  figure  6,  for  small  6  ,  a„„  can  be  considered 
equal  to  its  projection  on  the  1-3  plane  and  the  problem  is 
simplified  to  that  illustrated  in  figure  4,  where 

a33  =  90  '  a13' 
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The  cross  spectrum  between  wave  displacement  and  vertical 
velocity  for  this  situation  is  derived  in  Appendix  B  and 
can  be  written 


CCo(f)  -  i  Q(f))  df  = 


=  OlKw*  +  n  U  sina33)  ~  i(jl  W  cosa33) 


(32) 


Making  the  assumption  that  random  surface  displacements 
(e,g,  breaking  wave  crests)  are  essentially  uncorrelated 
with  turbulent  velocities  at  a  depth  of  5  meters  (the  minimum 


measurement  depth),  or  greater,  allows  the  n'w*  term  to  be 
neglected.   The  vertical  integral  scale,  measured  from 
temperature  data  to  test  this  assumption,  was  1,5  meters. 
The  cross  spectrum  is  now  defined  totally  by  wave-induced 
motions  and  the  phase  angle  in  the  wave  band  is 


Q(f)  df 


<{>  =  tan 


Co(f)  df 


n  W  cosa 
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n  U  sina 


(33) 


33 


Solving  for  a33: 


a33  =  tan 


-1 


W 


1  U  tan  <{> 


(34) 


43 


The  cross  spectra  between  waves  and  vertical  velocity  were 
obtained  for  the  two  shallowest  runs  and  cioo  was  calculated 
The  results  are  given  in  Table  III, 

Run  6:   seas  from  31QTj  swell  from  290T 
Run  10:  seas  from  30QT;  swell  from  270T 
vertical  flowmeter       Run  6  Run  10 

tilt  angle  Ca33)         3,54p  3,62° 


TABLE  III,   Tilt  angle  of  vertical  flowmeter  axis.   The 

azimuth  of  the  plane  of  vertical  and  horizontal 
flowmeter  components  was  270T, 


Using  a  value  of  a„„  equal  to  3,6  for  all  runs,  equation 
(18)  was  used  to  calculate  the  wave-induced  error  and  the 
correction  was  applied  to  the  measured  heat  flux.   The 
average  temperature  gradients  during  each  run  were  also 
measured  from  thermistors  approximately  one  half  meter  above 
and  below  the  point  of  flux  measurement.   The  eddy  diffusion 
coefficient  was  then  calculated  from 


6  'w' 
K„  =  Z-Z—  (35) 

U    36/8z 

The  turbulence  correlation  coefficient  (r*)  was  calculated 

©w 


from  corrected  values  of  0fw'  and  the  turbulent  intensities 
from  Tables  I  and  II.   The  results  are  given  in  Table  IV. 
The  turbulence  correlation  coefficients  were  all  less  than 
one  except  for  run  26.   A  turbulence  correlation  coefficient 
greater  than  one  indicates  that  there  were  errors  in  the 
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measurement  which  were  not  corrected  for  and,  hence,  run  26 
is  disregarded  in  the  discussion  below. 

The  variations  of  heat  flux  as  internal  waves  pass  a 
fixed  sensor  should  reflect  the  variations  in  velocity  and 
temperature  gradients  encountered.   In  the  absence  of  large 
scale  mixing  events,  such  as  the  breaking  of  internal  wave 
crests,  the  heat  flux  should  be  proportional  to  the  temper- 
ature and  velocity  gradients.   It  was  therefore  expected 
that  the  measured  heat  flux  would  increase  as  the  thermocline 
region  was  more  extensively  sampled.   Figure  7  and  Table  IV 
show  that  this  expectation  was  not  fully  realized,  since  a 
simple  relationship  between  temperature  gradient  and  heat 
flux  was  not  observed, 

Of  the  runs  conducted  at  night  (17,  18,  19),  the  first 

two  indicated  an  upward  transport  of  heat  presumably  due  to 

convective  overturn.   Both  these  runs,  although  close  to  the 

thermocline,  mainly  sampled  the  isothermal  regions  above  and 

below  it.   It  was  not  immediately  apparent  why  run  18,  below 

the  thermocline,  would  show  an  upward  heat  flux,  because 

surface  water  would  have  had  to  be  cooled  quite  radically  to 

sink  through  the  thermocline.   The  change  in  shape  of  the 

temperature  profile  before  and  after  the  run  indicated  that 

a  large  scale  mixing  event  might  have  occured  during  the  run, 

which  would  account  for  the  upward  heat  flux.   Run  19  sampled 

a  region  of  relatively  high  temperature  gradient  which 

resulted  in  a  net  downward  flux,  as  expected.   The  heat  flux 

was  less,  however,  than  for  a  daytime  run  in  a  region  of 

similar  temperature  gradient  (run  28a). 
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Runs  22  and  28  were  conducted  during  the  late  morning 
and  afternoon,  respectively.   Run  22  sampled  the  thermocline 
extensively  and  there  was  evidence  of  a  large  scale  mixing 
event,  both  in  the  temperature  time  series  and  in  the  BT 
samples.   The  heat  flux  was  consequently  large.   Run  28 
was  divided  into  two  runs  because  the  first  part  of  the  run 
sampled  an  isothermal  region  while  two  very  distinct  internal 
wave  troughs  passed  by  the  sensors  during  the  second  part  of 
the  run,   There  was  no  evidence  of  a  large  scale  mixing 
event  occuring  at  the  level  of  the  flux  measuring  sensors 
in  run  28,  although  the  thermistor  placed  a  meter  below  did 
record  a  possible  mixing  event.   This  run  was  considered  to 
be  measuring  only  small  scale  turbulent  flux  at  various 
distances  from  the  thermocline, 

The  high  turbulence  correlation  coefficients  observed  for 

runs  19,  22,  and  28a  correspond  to  relatively  long  periods 

of  time  during  which  the  sensors  were  in  the  thermocline 

region.   The  inference  can  be  made  that  gradient  type  eddy 

diffusion  in  this  region  is  an  efficient  heat  transport 

mechanism  which  can  probably  be  closely  approximated  by  a 

gradient  diffusion  model  of  the  form  of  equation  (35).   The 

relationship  between  K„  and  other  thermocline  parameters 

a 

(30/3z,  3U/8z)  could  not  be  determined  because  only  two  of 
the  six  runs  (19,  28a)  did  not  show  evidence  of  larger  scale 
mixing  events. 

On  a  qualitative  basis, the  results  showed  that  small 
scale  turbulent  heat  transport  in  the  upper  ocean  is  not 
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easily  modeled.    The  mechanism  for  transfer  seems  to  be  a 
combination  of  gradient  type  diffusion  and  convection.   Mixing 
at  large  scales  is  a  factor  and  may  be  due  to  breaking  of 
internal  wave  crests.   It  appears  that  large  scale  processes 
reestablish  sharp  temperature  gradients  which  are  again 
eroded  by  eddy  diffusion  and  convection.   Large  increases 
in  turbulent  kinetic  energy  were  not  observed  in  the  thermo- 
cline  region,  indicating  that  production  of  turbulence  by 
internal  wave-induced  shear  was  not  as  important  as  produc- 
tion at  the  surface.   The  local  Richardson  number  must  there- 
fore have  been  greater  than  the  critical  value  of  0.25  for 
the  maintenance  of  shear-induced  turbulence. 

VIII,   CONCLUSIONS 

Analysis  of  near  surface  velocity  records  in  shallow 
water  indicated  that  the  ratio  of  turbulent  to  wave-induced 
kinetic  energy  was  of  order  10"  .   Because  of  the  relatively 
large  wave-induced  fluctuations,  large  turbulent  flux 
measurement  errors  can  occur  for  small  errors  in  directional 
response,  phase  response,  or  instrument  alignment.   Alignment 
errors  are  by  far  the  most  troublesome  because  they  cannot 
be  determined  in  the  laboratory  before  measurements  are  made. 
Flowmeter  misalignment  of  3,6°  caused  errors  in  the  calculated 
heat  flux  of  about  25%  to  100%.   Momentum  flux  measurements 
would  have  been  affected  to  a  greater  degree.   It  was  shown 
how  momentum  flux  measurements  in  the  open  ocean  would  be 
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less  sensitive  to  misalignment  errors  than  shallow  water 
measurements.   Deep  water  misalignment  errors  would  be 
minimized  by  aligning  the  flowmeter  axes  as  near  as  possible 
to  the  predominant  direction  of  wave  motion, 

The  heat  flux  measurements  reported  here  were  only  an 
initial  effort  and  the  results  could  certainly  be  improved 
upon.   Verification  of  the  applicability  of  the  gradient  eddy 
diffusion  model  in  the  thermocline  was  not  possible  from  the 
sparse  data.   It  is  felt  that  analysis  of  the  remaining  NUC 
tower  runs  will  provide  sufficient  data  to  establish  a 
relationship  between  the  temperature  and  velocity  structure 
in  the  thermocline  and  the  eddy  diffusivity. 
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APPENDIX  A 

The  relative  magnitudes  of  the  wave-induced  and  turbulent 
fluctuations  can  be  approximated  by  assuming  that: 

1,)   the  wave-induced  and  turbulent  fluctuations  are 
independent ,    and 

2,)   the  wave-induced  fluctuations  can  be  described  by- 
linear  wave  theory, 

Using  the  first  assumption,  the  spectra  of  (1)  are  related 

according  to: 

Sff)  =  STT(f)  +  S  ,Cf)  (a) 


The  second  assumption  allows  S.T(_f)  to  be  written  in  terms  of 
the  spectrum  of  the  surface  displacement  (S  (f))  as: 

S0(f)  =  |H2(f)|  SnCf)  (b) 

where  H(f)  is  a  transfer  function  derived  from  linear  theory 
The  coherence  between  U  and  r\    is  given  by 


where  the  identity  holds  for  a  linear  constant  parameter 
system.   Now 

|SUn(f)|2     =     SyCf)     S^f)  (d) 

and  by  assumption  (1.) 

|STT  (f)|2  =  |S   (f)|2  =  STT(f)  S  (f)  (e) 
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The  measured  coherence  between  u  and  n   is 

2 

ylcf)  -i-aa. 


un 


suCf)  snCf) 


(f) 


substituting   (e)    and   Ca), 


2    ,_         SUCf)  Su(f)    ^   V(f) 


Y      (f)    = 

'un     ' 


s  Cf)  s  Cf) 

u  u 


(g) 


which  results   in 


S    .(f)    =   S    (f)    CI   -  Y      )    . 
u    „  u  un 
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APPENDIX  B 


1.   Cross  Spectrum  Between  Fourier  Components  of  r\   and  W* 


X?*=    ir^na-t   cot  o<       -f  Woo*  a*   owy 

-  £  W cm   fat  -o^J  +  i  Wccx  (cr-kf  0<3j) 

+  J:  7  U^vt^u+i  fat  +crr  -o<^) 


+  £  -  W  ^r.  z-i    004   (<-£   .  ^r   -PC33) 
*  a    7i   Kf-dUt  cr~t    Cot  (cr±  +o-t  -*-«*,3)J 
•  ci-b 


"*.  cz  n?  ~bv  co*  o<3j  ^on  C-crzj 


"4  J  W  <UH  ot^A^wt  cot  ATTj£r-L4l*G-T  M*ATTlr\\ 


^ 


Sfw*  Cp  v     O    />     J  *    *% 


OTHEAWt'sr.  . 
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2.   Cross  Spectrum  Between  Fourier  Components  of  U 


*   and  W* 


LT     —      CT *±i^t  o~t   com,,    +    Zlr  ooo  cr-t    cod«*3 


U*  = 


W> 


^Q.v,  (T)  -   TU-4  U'^UKfot+tx,,)^  (vt  +  vT+x^) 
*  4  U**u>^  Co-***,,)  M<n  (cr-k  ■+o-r-<Kz) 
•*    -jf  C/W^n  (v£  -kx,,  )  com  (<rt  +  art  -  <*„  ) 

-4       •£    JJlj^n  (<J±  +*,,)   CO^CrTt  +OZ  -r-«33) 

+  ■£  U  V«  Ccnf  -t<„)  <*m\  (o--i  rtrT  y.  oc^J 

*  if  U*  <d^n  (iri-tx,,  )  -don  fcr-t  +  o-zr  -<x,3  ) 
+  *£  UDJa^x  (at  -octt )  coc  (<3-£  *  err  -  <*s3) 

+  £  UK7  aJtx  (b~t  -  cxti  )  coo  (b~£  -a  erf  -^cx^3  ) 

•*•£  VJU  coo  (cri-Kx)   xi^t  (cr-t  -f-^t-f-^s) 
+  i  TaTCT  co<  Ccrii  -**,)  ^^  (*rt  -'err  -oc.j  ) 

•+  ~Y   W*   coo  (<r-t  -o<3l  J  coo  (cr*  v-o-r  -^3) 
•*■  ^  W  2  c<xj  fcrt  -<xr,/J  cos  (^  *<rr  +K33) 
-*  ^  £7  £7  c^j  fcrzi  ^,j  ^  Ccr^  *oT  ^a-m^ 
-*  -jj  TV  U  coo  (crt  -t-ocj,  )  sdtn  (vt  y-crr  -  <*,$) 

*  •£  TJ*  con  (a-±  -^  J  *tf  <Vii  +tr  T  -K^) 
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A    u    c°iC<xll+xls)eoi<rT    *  ^  H  ce&  (*u  -<*,3)  un  err 

•♦TV^COd    C«j(  +^33  )     +    +1*7*004    foe,- ©Cjj) 
CaJHEN  j?    -      ~T       '         O        OTH£*uJi's£ 

I 

SU'V*(p      -  4   V    C4KC*,,    C&SO<l2t     •*-   4    "W    CO"!  <Ksi    CO4  0<3i 

+  £-  I  4" UWcoitx,,  00^0^33    -  JUVC^I^,    C&KX^J 
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